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Abstract 

We show that the Lieb-Liniger model for one-dimensional bosons with repul- 
sive 5-function interaction can be rigorously derived via a scaling limit from a 
dilute three-dimensional Bose gas with arbitrary repulsive interaction potential 
of finite scattering length. For this purpose, we prove bounds on both the eigen- 
values and corresponding eigenfunctions of three-dimensional bosons in strongly 
elongated traps and relate them to the corresponding quantities in the Lieb- 
Liniger model. In particular, if both the scattering length a and the radius r of 
the cylindrical trap go to zero, the Lieb-Liniger model with coupling constant 
g ^ ajr^ is derived. Our bounds are uniform in g in the whole parameter range 
< g < oo, and apply to the Hamiltonian for three-dimensional bosons in a 
spectral window of size ^ above the ground state energy. 



1 Introduction 

Given the success of the Lieb-Liniger model [ISllII], both as a toy model in statistical 
mechanics and as a concrete model of dilute atomic gases in strongly elongated traps, 
it is worth investigating rigorously its connection to three-dimensional models with 
genuine particle interactions. A first step in this direction was taken in [17^ I16j . 
where it was shown that in an appropriate scaling limit the ground state energy of 
a dilute three-dimensional Bose gas is given by the ground state energy of the Lieb- 
Liniger modeL The purpose of this paper is to extend this result to excited energy 
eigenvalues and the corresponding eigenfunctions. 

The Lieb-Liniger model has recently received a lot of attention as a model for 
dilute Bose gases in strongly elongated traps [20lll2l[6l[9l[2l[l6l[2ll[23l[3]. Originally 
introduced as a toy model of a quantum many-body system, it has now become 
relevant for the description of actual quasi one-dimensional systems. The recent 
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advances in experimental techniques have made it possible to create such quasi one- 
dimensional systems in the laboratory [5l [18l [2H [25l [HI \7\. These provide a 
unique setting for studying matter under extreme conditions where quantum effects 
dominate. 

The Lieb-Liniger model describes n non-relativistic bosons in one spatial dimen- 
sion, interacting via a 5-function potential with strength g >0. In appropriate units, 
the Hamiltonian is given by 

n 

H^f' = Y.{-9f+^-'yHzi/i))+9 E ^(^-^j)- (1-1) 

i = l lS;*<i<" 

Here, we use the notation di = d/dzi for brevity. The trap is represented by the po- 
tential yll, which is assumed to be locally bounded and tend to infinity as \z\ oo. 
The scaling parameter £ is a measure of the size of the trap. Instead of the trap 
potential one can also confine the system to an interval of length I with appro- 
priate boundary conditions; in fact, periodic boundary conditions were considered in 

mm- 

The Lieb-Liniger Hamiltonian H^^'^ acts on totally symmetric wavefunctions 4> € 
L^(M"), i.e., square-integrable functions satisfying 4'{zi, . . . , Zn) = </>(z^(i), . . . , z^^^^-^) 
for any permutation vr. In the following, all wavefunctions will be considered sym- 
metric unless specified otherwise. 

In the case of periodic boundary conditions on the interval [0, £], Lieb and Liniger 
have shown that the spectrum and corresponding eigenfunctions of H^^ can be 
obtained via the Bethe ansatz [12j . In [Hj Lieb has specifically studied the excitation 
spectrum, which has an interesting two-branch structure. This structure has recently 
received a lot of attention [9l [3] in the physics literature. Our results below show 
that the excitation spectrum of the Lieb-Liniger model is a genuine property of dilute 
three-dimensional bosons in strongly elongated traps in an appropriate parameter 
regime. 

In the following, we shall consider dilute three-dimensional Bose gases in strongly 
elongated traps. Here, dilute means that ag^^^ <C 1, where a is the scattering length 
of the interaction potential, and g is the average particle density. Strongly elongated 
means that r <^ i, where r is the length scale of confinement in the directions 
perpendicular to z. We shall show that, for fixed n and i, the spectrum of a three 
dimensional Bose gas in an energy interval of size above the ground state 

energy is approximately equal to the spectrum of the Lieb-Liniger model (jl.ip as 
long as a <C r and r <^ i. The effective coupling parameter g in (II. ip is of the 
order g ~ a/r^ and can take any value in [0, cxo]. The same result applies to the 
corresponding eigenfunctions. They are approximately given by the corresponding 
eigenfunctions of the Lieb-Liniger Hamiltonian, multiplied by a product function of 
the variables orthogonal to z. The precise statement of our results will be given in 
the next section. 

We remark that the problem considered here is somewhat analogous to the 
one-dimensional behavior of atoms in extremely strong magnetic fields, where the 
Coulomb interaction behaves like an effective one-dimensional 5-potential when the 
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magnetic field shrinks the cyclotron radius of the electrons to zero. For such systems, 
the asymptotics of the ground state energy was studied in [Ij, and later the exci- 
tation spectrum and corresponding eigenfunctions were investigated in [3]. In this 
case, the effective one-dimensional potential can be obtained formally by integrat- 
ing out the variables transverse to the magnetic field in a suitably scaled Coulomb 
potential. Our case considered here is much more complicated, however. The cor- 
rect one-dimensional physics emerges only if the kinetic and potential parts of the 
Hamiltonian are considered together. 

2 Main Results 

Consider the Hamiltonian for n spinless bosons in three space dimensions, interacting 
via a pair-potential v. We shall write x = (x-'-,z) € M^, with x-*- G M^. Let 
G LJ^^(R^) and G LJ^j.(R) denote the (real- valued) confining potentials in x"*- 
in z-direction, respectively. Then 

n 

H^^^'''' = Y,{-^i + r-'V^{^i/r)+r'vHz,/£))+ J] a-\(|x, - x, |/a) . 
1=1 i<*<i<" 

(2.1) 

The trap potentials and confine the motion in the longitudinal (z) and the 
transversal (x^) directions, respectively, and are assumed to be locally bounded and 
tend to oo as \z\ and jx-*"! tend to oo. Without loss of generality, we can assume that 
yll > 0. The scaling parameters r and £ measure the size of the traps. 

The interaction potential v is assumed to be nonnegative, of finite range Rq and 
have scattering length 1; the scaled potential a~'^v{\ ■ \/a) then has scattering length 
a \14:\ ITSj and range aRo. We do not assume any smoothness or even integrability 
of V. In particular, we allow v to take the value +oo on a set of non-zero measure, 
corresponding to hard-sphere particles. 

Let and 6(x-'-) denote the ground state energy and the normalized ground 
state wave function of — A"*- -|- y-'-(x^), respectively. Note that 6 is a bounded and 
strictly positive function and, in particular, b G L^iM?) for any 2 < p < oo. Let 
also > denote the gap above the ground state energy of — + l/-'-(x^). The 
corresponding quantities for — A"*- -|- r~'^V-^{'x.^ /r) are then given by e^/r"^, jr^ 
and 6r(x^) = r~^6(x-'-/r), respectively. 

The eigenvalues of H^'l' ''^ will be denoted by E^^{n, r, £, a) , with A: = 1, 2, 3, ... . 
Moreover, the eigenvalues of H^f'^ in (jl.ip will be denoted by E^^ {n,£,g). Theorem 1 
shows that E^^{n,r,£,a) is approximately equal to E'[^{n,£, g) +e-'-/r^ for small a/r 
and r/i, for an appropriate value of the parameter g. In fact, g turns out to be given 

by 

g='-^[ |6(x^)|W. (2.2) 
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THEOREM 1. Let g be given in h2.S\) . There exist constants C > and D > 0, 
independent of a, r, i, n and k, such that the following bounds hold: 



(a) El^{n,r,£,a) > !^ + E'l^{n,i,g) {1 - n^) { I - ^EUn,i,g) ] (2.3) 



as long as E^^{n,(., g) < e^ jr^ , with 

'na\ 1/8 n /na\^/^ 



r]L = D 



2 C' 

V r 



(2.4) 



(b) El^{n,r,l,a)<— + E\^(nX9)i^-^^3r^ (2.5) 
whenever r/u < 1, where 



m-cD"\ (2,6) 



We note that E^^{n,£, g) is monotone increasing in g, and uniformly bounded 
in g for fixed k. In fact, Ef^{n,£, g) < Ef{n,£) for all g > and for all k, where 
Ef{n,£) are the eigenenergies of n non- interacting fermions in one dimension [8j. 
Using this property of uniform boundedness, we can obtain the following Corollary 
from Theorem [T] 

COROLLARY 1. Fix k, n and i. If r ^ 0, a ^ in such a way that a/r — > 0, 
then 

,^EMhl^A.ZJ^ = l, (2.7) 

El{n,l,g) ^ ' 

Note that by simple scaling E^^{n,r,£,a) = £~'^E^^{n,r/£,l,a/i), and likewise 
E'^^{n,£, g) = £~'^E^^{n,l,£g). Hence it is no restriction to fix £ in the limit con- 
sidered in Corollary [H In fact, £ could be set equal to 1 without loss of generality. 
Note also that the convergence stated in Corollary [1] is uniform in g, in the sense 
that g ~ ar~^ is allowed to go to +00 as r ^ 0, a ^ 0, as long as gr ~ a/r — > 0. 

In order to state our results on the corresponding eigenfunctions of H^^'^'"" and 
^id '^ ' we first have to introduce some additional notation to take into account the 
possible degeneracies of the eigenvalues. Let ki = 1, and let ki be recursively defined 

by 

ki = mm[k : E'[^in,£,g) > E';^-'{n,£,g)] . 

Then E'l^{n,£,g) < E\'+\n,£,g), while ^(n, 5) = i?fA(n,^,g) for < j < 
fej+i — ki. That is, k counts the energy levels including multiplicities, while i counts 
the levels without multiplicities. Hence, if ki < k < ki^i, the energy eigenvalue 
E'^^{n,£, g) is /cj+i — ki fold degenerate. Note that ki depends on n, £ and g, of 
course, but we suppress this dependence in the notation for simplicity. 
Our main result concerning the eigenfunctions of H^^'^'"" is as follows. 
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THEOREM 2. Let g be given in 112. Let be an eigenfunction of H^^'^'"" with 
eigenvalue E^^{n,r,i,a), and let tjji (ki < I < /cj+ij be orthonormal eigenf unctions 
of H^f'^ corresponding to eigenvalue Ef^{n,£,g). Then 



ki+i-l 

E 



i=l 



-1 



> 1 



EtT'E{^in,e,g) 



E':'+'{n,e,g)-E'Un,e,g) El:,{n,£,g) - Ely\n,£,g) 



+ 



{rye^)E>l^{n,£,g)) 
(1 -?7u)(l -r/L) 

E-U'El^in,£,g) 



ki-l. 



as long as r/u < 1, ??l < 1 CLnd E'[^{n,£,g) < e^/r^. 

In particular, this shows that ^'^(xi, . . . , x„) is approximately of the product 
form ipk{zi, . . . , Zn) Yii ^r(x^) for small r/£ and a/r, where ipk is an eigenfunction of 



H^^'^ with eigenvalue E^^{n,£,g). Note that although ^'^ is close to such a product 
in L^(M^"') sense, it is certainly not close to a product in a stronger norm involving 
the energy. In fact, a product wavefunction will have infinite energy if the interaction 
potential v contains a hard core; in any case, its energy will be too big and not related 
to the scattering length of v at all. 



COROLLARY 2. Fix k, n, t and go > 

eigenspace of H^^'^^ with eigenvalue Ef^{n,£, go] 
onto the function IlILi^K^) ^ If r 

g = 87r|[6|[4a/r-^ go, then 



Let Pg^ denote the projection onto the 

, and let denote the projection 
^ 0, a — > m such a way that 



lim(^'JPi 



1. 



(2.9) 



p2n\ 



Here, the tensor product refers to the decomposition L"^ 
into the transversal (x^) and longitudinal (z) variables. 

We note that Corollary [2] holds also in case go = co. 
has to be defined as the spectral projection with respect to the limiting energies 
\mig^rx> E^^{n,£, g). Using compactness, it is in fact easy to see that the limit of 
Pg^ Id as (7o — > oo exists in the operator norm topology. We omit the details. 

Our results in particular imply that Hl^^'^'"" converges to H^^'^ in a certain norm- 
resolvent sense. In fact, if a — > and r — > in such a way that g = 87r||6|j|a/r^ —>■ go, 
and A G C \ [El^{n,£, go), oo) is fixed, then it follows easily from Corollaries [1] and [2] 
that 



In this case, P^^ 



lim 



1 



A + ne^/r'^ 



H. 



n,r,i,a 
3d 



1 



X-H- 



Id 



0. 



Our main results. Theorems [T] and O can be extended in several ways, as we will 
explain now. For simplicity and transparency, we shall not formulate the proofs in 
the most general setting. 
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• Instead of allowing the whole of as the configuration space for the vari- 
ables, one could restrict it to a subset, with appropriate boundary conditions 
for the Laplacian A"*-. For instance, if V'^ is zero and the motion is restricted 
to a disk with Dirichlet boundary conditions, the corresponding ground state 
function b in the transversal directions is given by a Bessel function. 

• Similarly, instead of taking M as the configuration space for the z variables, one 
can work on an interval with appropriate boundary conditions. In particular, 
the case V^^ = with periodic boundary conditions on [0, i] can be considered, 
which is the special case studied by Lieb and Liniger in [12^ 111]. 

• As noted in previous works on dilute Bose gases [151 [T^ , the restriction of v 
having a finite range can be dropped. Corollaries [U and [2] remain true for all 
repulsive interaction potentials v with finite scattering length. Also Theorems[T] 
and [2] remain valid, with possibly modified error terms, however, depending on 
the rate of decay of v at infinity. We refer to [H] for details. 

• Our results can be extended to any symmetry type of the ID wavefunctions, 
not just symmetric ones. In particular, one can allow the particles to have 
internal degrees of freedom, like spin. 

• As mentioned in the Introduction, in the special case k = 1, i.e., for the ground 
state energy, similar bounds as in Theorem 1 have been obtained in |T^. In 
spite of the fact that the error terms are not uniform in the particle number 
n, these bounds have then been used to estimate the ground state energy in 
the thermodynamic limit, using the technique of Dirichlet-Neumann bracket- 
ing. Combining this technique with the results of Theorem 1, one can obtain 
appropriate bounds on the free energy at positive temperature and other ther- 
modynamic potentials in the thermodynamic limit as well. In fact, since our 
energy bounds apply to all (low- lying) energy eigenvalues, bounds on the free 
energy in a finite volume are readily obtained. The technique employed in [T7] 
then allows for an extension of these bounds to infinite volume (at fixed particle 
density) . 

In the following, we shall give the proof of Theorems [T] and [2j The next Section [3] 
gives the proof of the upper bound to the energies E^^{n,r,£,a), as stated in Theo- 
rem [T]^b). The corresponding lower bounds in Theorem [T]^a) are proved in Section HI 
Finally, the proof of Theorem [2] will be given in Section [5j 

3 Upper Bounds 

This section contains the proof of the upper bounds to the 3D energies E^^{n, r, £, a), 
stated in Theorem [2b). Our strategy is similar to the one in [n\ Sect. 3.1] and we 
will use some of the estimates derived there. The main improvements presented here 
concern the extension to excited energy eigenvalues, and the derivation of a bound 
that is uniform in the effective coupling constant g in (|2.2p . In contrast, the upper 
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bound in [T71 Thm. 3.1] applies only to the ground state energy E^^{n,r,i,a), and 
is not uniform in g for large g. 

Before proving the upper bounds to the 3D energies E^^{n,r,£,a), we will prove 
a simple lemma that will turn out to be useful in the following. 

LEMMA 1. Let H be a non-negative Hamiltonian on a Hilbertspace 7i, with eigen- 
values < El < E2 < E^ < . . . . For k > 1, let fi, . . . , fk Ti.. If, for any {aj} with 
Yli=i l^iP — 1 '^^ have 



k ^ 
J2i=iO'ifi 



> 1 — e for e < 1 



and 

{Ei=lO'if^\H\E^=la^f^)<E, 

then Ek < E{1 - e)-i. 

Proof. For any / in the fc-dimensional subspace spanned by the /j, we have 

{f\H\f)<E{l-e)-\f\f) 

by assumption. Hence E/^. < E/{1 — e) by the variational principle. ■ 

Pick {oj} with Yli=i — 1' ^'^d let = Yli=i ^i''Pi^ where the tf^i are orthonor- 
mal eigenfunctions of H^^ with eigenvalues E^^^{n,£, g). Consider a 3D trial wave 
function of the form 

n 

^>(X1, . . . ,X„) = 0(Z1, . . . , Z„)F(X1, . . . , X„) 6r(Xfc ) , 

k=l 

where F is defined by 



F(xi,...,x„) = JJ/([xi - Xjl) . 



i<j 

Here, / is a function with < / < 1, monotone increasing, such that f{t) = 1 for 
i > i? for some R > aRo. For t < R we shall choose f{t) = fo{t)/ fo{R), where /o is 
the solution to the zero-energy scattering equation for a~'^v{\ ■ \/a) [T^llSj. That is, 

normalized such that lim(^oo fo{t) = 1- This /o has the properties that fo{t) = 1—a/t 
for t > aRo, and /o(i) < min{l, a/t}. 

To be able to apply Lemma [H we need a lower bound on the norm of ^. This can 

f 2") 

be obtained in the same way as in [17, Eq. (3.9)]. Let denote the two-particle 
density of (j), normalized as 

J2) 



{z, z')dzdz' 
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Since F is 1 if no pair of particles is closer together than a distance R, we can estimate 
the norm of ^ by 

n{n — 1) 



(^>|$) > 1 



> 1 



n(n — 1) 



(3.1) 



2 

where Young's inequality [13^ Thm. 4.2] has been used in the last step. Since F < 1, 
the norm of ^ is less than 1, and hence 



1 > ($1$) > 



r 



(3.2) 



Next, we will derive an upper bound on {^\H.^^ ' ' |<I>). Define G by <I> = GF. 
Using partial integration and the fact that F is real-valued, we have 

(Cl>|-A,|cl>) = - 



F'GAjG+ / \G\'\VjF\\ 



Using —A-^br = {e^ /r'^)br — r '^V-^{ ■ /r)bj., we therefore get 



($|//3"/'^'"|ci>) = I^($|cI>) + 



F'\{br{^fn{H\ 



+ / \G? (ElV,i^l' + E«"V|x,,-x,I/a)|Fp 



(3.3) 



With the aid of Schwarz's inequality for the integration over the z variables, as 
well as F < 1, 



1/2 



1/2 



(3.4) 



Here, we have used that </> is a linear combination of the first k eigenfunctions of 
H^^'^ to obtain the final inequality. The term in the second line of (j3.3p is bounded 

by 



n{n — 1) 



.(2) 



(z, z)dz ( 1 



n{n — 1) ttR"^ 



14 ; ) 
(3.5) 
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as an argument similar to (j3.ip shows. 

The remaining last term in (j3.3p can be bounded in the similar way as in \17\ 
Sect. 3]. In fact, by repeating the analysis in [171 Eqs. (3.12)-(3.19)] 



< 



n{n — 1) 



.(2) 



{z, z')h{z — z')dzdz' 



2 , -Jl^ll^ ll^lll 



+ -n(n - l)(n - 2)- 



.(2) 



(z, z')m{z — z')dzdz' . (3.6) 



Here, 



and 



/.(z) = / (/'(|x|)2 + ia-\(|x|/a)/(|x|)2) d^x 



m(z) = / /'(IxDd^xJ 



(The function m was called A; in [17].) Note that both /i and m are supported in 
[-R,R], and J^h{z)dz = 47ra(l - a/R)~^. 

We now proceed differently than in Our analysis here has the advantage of 
yielding an upper bound that is uniform in g. Let ip G H^(W). Then 



'm2| 



dt 



dt 



<2\^iz')\ 



difit) 



dt 



<2\y.{z')\\z-z'\'/' 



dt + 2 



dip 



dip{t) 



dt 



dt 



(3.7) 



dz 



1/2 



+ 2\z - z'\ 



dip 



dz 



Here, we used \(p{t)\ < \ipiz') \ + f^, ^^^^ dt for z' < t < z for the first inequality and 
applied Schwarz's inequality for the second. 

We apply the bound (13. 7p to (z,z') for fixed z'. Using the fact that the 
support of h is contained in [—R, R], we get 



g\ {z, z')h{z — z')dzdz' — / h{z)dz I q\ {z,z)dz 

Jr Jr 

<2i?V2 I h{z)dz j j d,^Qf{z,z') dz\ ^ef{z',z')dz' 

+ 2R f h{z)dz [ dJgfiz.z') 
Jr Jr2 ^ ^ 



<-EUn,l,g) / h{z)dz 



2R 



dzdz 

1/2 



(n - l)g 



+ R 
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('2) 

Here, we used Schwarz's inequality and the fact that g^^nin — 1) j^Q^^ {z,z)dz < 
E\^{n4,g) and 



dzdz' < ( 



dzf 



(t>) <EUnA.9)ln. 



The same argument is used with h replaced by m. Now h{z)dz = 47ra(l 
a/R)-^, and ^IZt Eq. (3.22)] 



Imlino < 



2'Ka 



1 - a/R 



(1 + ln(i?/a)) 



[ 2TTaR ( a 
/ m{z)dz < I 1 



1 - ajR \ 2R 



Therefore 



_ ra(n-l) l + K 

IM < ^ ^ ' g- 



2 ''I -a/R 
where we denoted 



il\z^z)dz + - 



n 



2R 



[n-l)g 



1/2 



+ R 



Efdin,i,g) 



(3.8) 



27r, n l + ln(R/a) ,,,,,2 

Putting together the bounds p.4p . (|3.5p and p.Sp . and using again the fact that 
g^n{n — 1) Q^^\zi z)dz < E^^{n,i, g), we obtain the upper bound 



h: 



n,r,£,a ^6 
3d Z2 



+ 



1 + iC 
1 - a/R 



(1 - a/R) 
{2Rg{n-l)Y/^ + {n-l)Rg 



It remains to choose R. If we choose 



R-" 



ar 



then 



^3d 



ne 



<£:fd(n,^,ff) 1 + C 



for some constant C > 0. Moreover, from ()3.2p we see that 

m>i-c(!^)''\ 

Hence the upper bound (j2.5p of Theorem [1] follows with the aid of Lemma [TJ 
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4 Lower Bounds 



In this section, we will derive a lower bound on the operator Hnf'^'"" in terms of 



the Lieb-Liniger Hamiltonian H^^'^. In particular, this will prove the desired lower 
bounds on the energies E^^{n,r,i,a). Our method is based on [17, Thm. 3.1], but 
extends it in several important ways which we shall explain. 

Let $ be a normalized wavefunction in We define / G by 

„ n 
f{zi,...,Zn) = $(xi,...,X„)rr6r(x^)d^X^. (4.1) 

Moreover, we define F by 

n 

$(xi, . . . ,x„) = F(xi, . . . ,x„) JJ 6,.(x^) . 



fc=i 



Note that F is well-defined, since br is a strictly positive function. Finally, let G be 
given by 



G(xi,. . . ,X„) = ^>(xi,. . . ,X„) - /(Zi,. . . ,Z„) 6r(Xfc ) . 



(4.2) 



k=l 



Using partial integration and the eigenvalue equation for br, we obtain 



H. 



ne 



3d 



i=l 



ViF|2+rVll(zi/^)|F|2 (4.3) 

n 

n^r(x^)2d3xfc. 



+ i ^ a ^v{\-Ki-^j\/a)\F\- 

Now choose some R > aRo, and let 

3{R^ - aR^)-^ for aRo < r < R 



k=l 







otherwise 



U{r) = 

For (5 > define Bs C by 

Bs = |x^ G : 6(x^)2 > ^1 . 

Proceeding along the same lines as in [171 Eqs. (3.31)-(3.36)], we conclude that, for 
any < e < 1, 



n 

j=i 
+ 



e|V,^F|^ + a'C/(|x,-x,(,)|)xB,(x^(,)/r)|F| 



k=l 



eldi^l"^ + (1 - e)\d^<l>\'^Xrnink |x,-x,|>i?(xi) 



n ^'^fc 



(4.4) 



k=l 
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Here, x^^j) denote the nearest neighbor of Xj among the Xj with j ^ 
_L-component, and a' = a(l — — 2i?||V6^||oo/r(5). (It is easy to see that is 
a bounded function; see, e.g., the proof of Lemma 1 in the Appendix in [IT].) The 
characteristic function of Bs is denoted by XBg^ ™d Xminfe [x,-Xfe[>_R restricts the Xj 
integration to the complement of the balls of radius R centered at the x^ for k ^ i. 
That is, for the lower bound in (j4.4p only the kinetic energy inside these balls gets 
used. (Compared with [TTl Eq. (3.36)], we have dropped part of the kinetic energy 
in the x^ direction in the last line, which is legitimate for a lower bound.) 

For a lower bound, the characteristic function Xmink\x^-xk\>R could be replaced 
by the smaller quantity Xmiiife |zj-Zfc|>i?; ^ was done in [17]. We do not do this here, 
however, and this point will be important in the following. In particular, it allows 
us to have the full kinetic energy (in the z direction) at our disposal in the effective 
one-dimensional problem that is obtained after integrating out the x^ variables. In 
contrast, only the kinetic energy in the regions \zi — Zk\ > R was used in [17j to 
derive a lower bound on the ground state energy. The improved method presented 
in the following leads to an operator lower bound, however. 

We now give a lower bound on the two terms on the right side of (j4.4p . We start 
with the second term, which is bounded from below by 

„ n 

„ n 
- (1 - e) / |ai$Pxniin, \z,-z^\<R{zi)Xmin^ |x^ ~x^ | <fl(x»^ ) H ^^""^ ' 

To estimate the last term from below, consider first the integral over the x^ variables 
for fixed values of zi , . . . , 2:„ . Using (j4.2p , this integral equals 

r 2 " 

/ 5i/nLl^(Xfc ) + diG Xmin, \^i-^j:\<R{^i) H ^^^^ 

JR2" ' k=l 

« n 
< (1 + r?-l) / Xrn^n, |x^-x^|<ii(x.^) H ^ri^ifcP^i 

„ n 

+ (1 + 7?) / |a,G|2J]d2^^ 

for any 77 > 0, by Schwarz's inequality. It is easy to see that 

nin — 1) ttR'^ 



— Il"ll4 

k=l 



using Young's inequality, as in (jS.ip . 
Now 



n „ ra „ n 

fc=i ^R^" t=i ^K^" fc=i 
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since, by the definition of G in (j4.2p . diG is ortliogonal to any function of the form 
^(2:1, . . . , Zn) Y\k ^r(x^). We have thus shown the the second term on the right side 
of (j4.4p satisfies the lower bound 



/ 



e\dM^ + (1 - e)|a,«>pxmmfc |x.-x,|>fl(xi) 



k=l 



/ / s / i\ n(n — 1) TiE? „, „A 



k=l 



+ {l-(l-e){l + rj)) [ \diG\'f\d^^k 



(4.5) 



for any 77 > 0. We are going to choose rj small enough such that (1 — e) (1 + r/) < 1, in 
which case the last term is non-negative and can be dropped for a lower bound. Note 
that the right side of ()4.5p contains the full kinetic energy of the function /. Had 
we replaced Xminfe |xi-Xfe|>i? by the smaller quantity Xmiuk \z,-Zk\>R the left side of 
(|4.5p . as in [17], only the kinetic energy for jzj — Zfc| > i? would be at our disposal. 
This is sufficient for a lower bound on the ground state energy, as in [17], but would 
not lead to the desired operator lower bound which is derived in this section. 

We now investigate the first term on the right side of (j4.4p . Consider, for fixed 
zi, . . . , Zn, the expression 



n „ n 

Proceeding as in Eqs. (3.39)-(3.46)], we conclude that 



(4.6) 



„ n 

(HD >a"Vd(z, -z,) / |$(xi,...,x„)|2TT(i2x: 



k ; 



where 



a' 1 



3n ar 



ee^ 1 - {aRo/R)^ [ ee^ R 



1 



1 - {aRo/Rf 



X (l-(n-2)^||6||L 



(4.7) 



(which is called a'" in [I?]), and 

d{z -z')= [ 6,(x^)26,(y^)2[/(|x - y|)xe.(y^A)d'x^(i2y^ • 

Note that d{z) = if |z| > R. Recall that the quantity e"*" in (j4.7p is the gap above 
the ground state energy of — A"*- + V^, which is strictly positive. 
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By using 

„ n 



k=l 



Let g' = a" d{z)dz. It was shown in [T71 Eqs. (3.49)-(3.53)] that, for any x > 
and if e H'^ l^), 



a" I d{z — z')\lp{z)'^ dz 



> g' max MzM^ 1 



2g'R 



K / \dip'\^dz 



for fixed z' . In particular, applying this estimate to / for fixed Zj, j 7^ i, 

„ n 
a"V / d{zi-Zj)\j\^T\dzk 

„ n n „ n 

>g"y2 s{zi-z,)\f\^Y\dzk-{n-i)^y2 i^^/i'TT^^'^^ 

T^i^^" k=i ^-^iK" fc=i 



with 



Hence we have shown that 



2g'R 



(4,8) 



E 



e|V,^F|^ + a'[/(lx,-Xfc(,)|)x6,(x^(,)/r)|F 



„ n n „ n 

>/V/ 6{zi-z,)\f\^T\dzk-{n-l)>,y^ mi'Yldzk. (4.9) 

This concludes our bound on the first term on the right side of (j4.4p . 

For the term involving V^^ in (j4.3p . we can again use the orthogonality properties 
of G in (j4.2p . as well as the assumed positivity of V\ to conclude that 

n „ n 

V / rVii(z,/^)|ci>(xi,...,x„)pn^'^'= 

n „ n 
>V/ r'vHzi/i)\f{zi,...,Zn)\'T\dZk. (4.10) 

By combining (fO|) and (fO|) with the estimates (fOj) and (|4.10p . we thus 
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obtain 



h: 



n,r,l,a 



3d 



> 1 - (1 - e) (1 + 7?-!) - (n 



+ 



1=1 



i<3 



[ \dif\^f]dZk 

n 



(4.11) 



k=l 



Recall that g" is given in (j4.8p . This estimate is valid for allO<e< 1, x>0 and 
r] > such that (1 — e)(l + ry) < 1 in order to be able to drop the last term in (|4.5p . 

To complete the estimate, it remains to give a lower bound to J^d{z)dz. As in 
[HI Eq. (3.48)], we can use \b{x-^f - b{y^)'^\ < R\\\/b'^\\oo for jx-^ - y-^\ < R to 
estimate 



/ d{z)dz >%( [ b{^^fd^^^ - i?||V62 



47r 



This leads to the bound 



> 1 



2gR 



> — (l|fe|l!-'^-^l|V6'||oo/r) 



R\\Vb'^\\oo/r 



[1 - e) (1 - 2R\\Vb^\\oo/r6) 1 - (n - 2) 



ttR' 



M 

2 ll"lloo 



X 1 



3n ar 



1 



ee^R " "^1 - (ai?o/^)^ 



It remains to choose the free parameters R, 6, x, e and r/. If we choose 



e = T] 



1 /na\5/i2 
n V r 



Eq. (|4.1ip implies 



if: 



ne 



3d 



r / 



/ 



for some constant Z) > 0. In particular, if denotes the projection onto the 
function 0^=1 ^r'(x^) G L^(M^"), we have the operator inequality 



-^3d - — > (1 - Vl) Pr ®H^' 



'Id 



in the sense of quadratic forms. (Here, r^L is defined in ([2 
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Recall that that the gap above the ground state energy 
is e-^ /r"^. Hence, using the positivity of both V^^ and 

In particular, 

^sT'''^ - ^ > (1 - 7) (1 - r/L) Pr^ ® ^r/'' + 7^ (l - Pr"-) ® 1 (4.12) 

for any < 7 < 1. If we choose 7 = (r^/e-'-)£'^^(n,^, 5), the lowest k eigenvalues 
of the operator on the right side of (|4.12p are given by (1 — 7)(1 — r]i^)El^{n,£, g), 
1 < i < k. This implies the lower bound (j2.3p of Theorem [TJ 



5 Estimates on Eigenfunctions 



In this final section, we shall show how the estimates in the previous two sections 

Tn,e,g 



can be combined to yield a relation between the eigenfunctions 4' a,- of Hnf' and 



the eigenfunctions il^k of -f^ij • We start with the following simple Lemma. 

LEMMA 2. Let H be a non-negative Hamiltonian on a Hilbertspace Ti, with eigen- 
values < E'l < E2 < ... and corresponding (orthonormal) eigenstates tpi. For 
k > 1, let fi, 1 <: i < k, be orthonormal states in 7i, with the property that 
{fi\H\fi) < r]Ei for some rj > 1. If E^+i > Ej,, then 

. , J^k+l — 

Proof. Let P^ be the projection onto the first k eigenfunctions of H. Then H > 
HPk + Ek+i{l — Pk)- By the variational principle, 

k k 

Y,{fi\HPk + Ek{l - Pk)\fi) >Y.Ei, 

i=l i=l 

and hence 

k k k 

1=1 i=l i=l 

By assumption, {fi\H\fi) < rjEi, implying the desired bound 



i=l 
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COROLLARY 3. Assume in addition that Ei^i > Ei for some 1 < I < k. Then 



(ri - 1) EU E^ (r? - 1) E:=i 



i,j=l+i 



Ek+i — Ek 



Ei+i - El 



Proof. Applying Lemma [2] and using the orthonormality of the ipi and the fi, we 
have 

k k I 

i,j=l+l i,j=l i,j=l 

Ik k I 

i=l j=l i=l j=l 

>,_ c»:i)Etig. ^,_ (')-i)sLi^. _,,, 



Ek+i — Ek 



Ei+i — El 



Let again denote the projection onto the function Y[k=i^r{^k) ^ L'^{M?"') 
and recall inequality (I4.12p . which states that 



Hsf'^ - ^ > (1 - 7) (1 - Vl) + 7^ (l - Pr^) ® 1 (5-1) 

for any < 7 < 1. As already noted, the choice 

^ = {rye^)EUnJ,g) 

implies that the lowest k eigenvalues of the operator on the right side of (|5.1|) are 
given by (1 — 7)(1 — r]i^)El^{n,£, g), 1 < i < k. Here, we have to assume that 
~±^j,2 y E^^{n,l,g) in order for 7 not to be greater than one. The corresponding 
eigenfunctions are 



'4)i{zi,...,Zn)WK{y.i 



k=l 



We can now apply Corollary [31 with H equal to the operator on the right side of 
dSlD, fi = ^i,l = fei-l and k = k,+i-l. Since El^{n,r,i,a) < E'l^{n,i,g)/{l-r]v) 
by Theorem 1(b), we conclude that 



i=l 



> k. 



i+i 



ki 



r\^k . 



1 



(l-7)(l-7?u)(l-r?L) 



E- 



Id 



E' 



+ 



Id 



^id - ^id 
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as long as 7 < 1, ryu < 1 and ?7l < 1- Here, we used the notation E^^ = El^{n, i, g) 
for short. Since 



for every k, this imphes Theorem [21 
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